





LASER PULSE ANALYSIS 


C. G. Lehr, C. R- H. Tsiang, G. M. Mendes, and R. J. Eidred 


/ 




To be presented at the First International Symposium 
on Use of Artificial Satellites 
for Geodesy and Geodynamics 


Athens, 14-21 May 1973 


N75-17657 


ft) ASA-CR-1 36759) LA SEE PULSE ANALYSIS 
(Smithsonian A strophysical observatory)^ ^ 

14 p HC $3.25 Unclas 

_ , G 3/3.6 11651 



Smithsonian Institution 
Astrophysical Observatory 
Cambridge, Massachusetts 02138 

■j) **.l V- W 


305-037 



LASER PULSE ANALYSIS 


C. G. Lehr, C. R. H. Tsiang, G. M. Mendes, and R. J. Eldred 


1. INTRODUCTION 

In a pulsed laser ranging system, the range error for a single-electron return is 
equal to the standard deviation of a probability distribution having the shape of the 
transmitted pulse. A typical Q-switched ruby laser transmits a pulse whose shape is 
closely appi'oximated by a normal-density distribution 

= a _1 ( 2ir)~ 1//2 exp (-t 2 /2 a 2 ) , (1) 

where the center of the pulse is at t = 0 and the pulse's standard deviation is a. For 

a typical pulse, the value of t, the time between half-intensity points, is about 20 ns, 

1/2 

and the corresponding cr value is t/ 2 (2in 2) ' or cr = 8. 5 ns. When a single photo- 
electron is received, the range error is this cr value times half the velocity of light, 

0. 15 m/ns, or 1. 3 m. If, instead of a single electron, the return contains N electrons, 
the returned pulse gives N independent measurements of the range. The average of 

these measurements has a range error that is smaller than the single-electron meas- 

1/2 

urement by a factor of N 7 . Hence, it can be seen that the range error of a 20-ns laser 
can be less than 0. 5 m when there are 7 or more electrons in the return and less than 
0.1m when the number of electrons is 163 or more. 

It is seldom practical to resolve each electron in the return, determine its respec- 
tive range, and compute the average. Instead, the digital range counter stops at a point 
on the pulse's leading edge, this point being determined by the threshold- voltage setting 
of the input discriminator of the counter's STOP circuit. A series of incremental 
ranges referenced from this threshold point are then measured to various points on the 
pulse, and these points are referred to the center of the pulse. If there is a significant 
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r^je-to-pulse variation in the shape or size of the transmitted pulse, a similar cor- 
rection for the time between the START threshold and the transmitted pulse center can 
be obtained by the same procedure. 

At the observing stations of the Smithsonian Astrophysical Observatoiy (SAO), 
oscilloscopic presentations of the transmitted and received pulses have been photo- 
graphed in a study of the laser system's ranging accuracy. Figure 1 is a block 
diagram of the range-measuring equipment. When the counter's START and STOP 
thresholds are reached, the sweep circuits of oscilloscope channels 1 and 2, respectively, 
are triggered, and the transmitted and received pulses are displayed and photographed. 

Figure 2 shows typical transmitted and received pulses on a time axis, with the 
amplitude of the transmitted pulse greatly attenuated. The counter measures the time 
T c between the threshold points T TT and T RT - Since the corrected range is the time 
T R between the pulse centers T TC and T RC , a correction A must be added to T c , 
where 

A = ( t rc " T RT^ “ * T TC " T TT^ ‘ ^ 

This correction is known once the threshold points and pulse centers have been deter- 
mined. Figure 3 shows a typical pair of transmitted and received pulses as displayed 
on an oscilloscope photograph, measured, and plotted by a digital computer. 

This paper is concerned with methods for locating the points T,^, T^, T R ,p, 
and on presentations like those of Figure 3. There are errors involved in the 
determination of each these quantities, and an attempt is made to separate their effects 
on the overall range correction. Several series of corrected range measurements for 
fixed reflectors and satellites were obtained. Residuals were computed by fitting the 
range measurements to either fixed- reflector distances or short arcs of satellite 
orbits. Rms values of these residuals are presented but not enough data have yet been 
analyzed to be sure that the results are either typical or optimum. 
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2. THE EXPERIMENTAL DATA 


The oscilloscope photographs of the transmitted and received pulses are obtained 
with a Tektronix model 454 oscilloscope, whose bandwidth is 150 MHz. The pulses 
were originally photographed on Polaroid film; later, 35-mm cameras were used. At 
first, the x-axis sweep was 50 ns, but this was later changed to 100 ns to avoid the 
truncation of the leading and trailing edges of the pulses. The Polaroid films were 
measured by hand at 50 points uniformly spaced along the x axis. The 35-mm films 
were measured at 100 time points by projecting the film on a screen and using a meas- 
uring machine. The measurements were punched on cards, read into a digital computer, 

* 

and stored on magnetic tape. Two sets of the data used in this report were obtained 
from only 30 to 40 measured points per pulse. These were converted to 100 data points 
by constructing a complete Fourier series of 15 harmonics. Tests showed that data so 
obtained are nearly identical to those obtained by the other method. The data for each 
frame are stored on magnetic tape as 200 amplitudes for 100 equally spaced time points 
on the transmitted and received pulses. For the 50-ns sweep, the points were 0. 5 ns 
apart, whereas for the 100-ns sweep, they were 1 ns apart. 

3. ANALYSIS OF THE EXPERIMENTAL DATA 

In computing the range correction A, the transmitted and received pulses must be 
located in time with respect to the starting and the stopping of the range counter. 

These locations must be determined to within a small fraction of the pulse durations. 
Hence, some point on the received pulse must be identified with a similar point on the 
transmitted pulse. The nearly normal shape of the transmitted pulse makes its "center" 
an obvious choice for this point. However, as will be shown below, the corresponding 
center of the irregular received pulse can be specified in several different ways. It is 
also possible to align the 100-point sets of the transmitted and received pulses by using 
a correlation technique that does not require the pulse center to be identified. 

One way to find the center of the pulse is to compute its centroid. If the pulse 
shape defines a density distribution on the time axis, then the centroid is the mean 


The fixed-reflector returns at the station in Arizona and the BE-C satellite returns 
at the station in Peru. 
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value of this distribution. If RP(t) is a continuous function representing the received 
pulse and t is its centroid, then the probability density distribution is 


oo 

f(t) * RP(t )/j RP(t) 

-OO 


dt 


and 


t C *J tf(t)dt 

-oo 


( 3 ) 


When the integral is evaluated on the computer, summation, rather than integration, 
is used. 

The centroid determination does not explicitly use the fact that the transmitted 
pulse and the probability- density function producing the received pulse are symmetric 
about their mean values. To find this symmetric center for an irregularly shaped 
received pulse, we set its time origin (t' = 0) at successive points t^tg, . . . ,t., . . . ,t i0Q . 
For each of these origins, we write fit) as a sum of an even and an odd function, viz. 


fito = f e (t') + f 0 (t') , (4) 

where f 0 {t') and f Q (t') can be determined as follows: 

f e (t') = | [fit') + f(-t')] , (5a) 

f Q (t') = | [fit') - f(-t')] . (5b) 


The symmetric center, t g , is defined as that value of t. that minimizes the odd func- 
tion in the sense that 


Here we are neglecting the distortion that is due to the nonsymmetric configurations 
of retroreflectors on satellites currently in orbit. Since the amount of this distortion 
varies as the satellite transits, it cannot readily be taken into account by the methods 
we are now using. 
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( 6 ) 


J l f 0 (t')! dt'= min . 

0 

In practice, the t. values are taken at 0. 5- or 1-ns intervals (depending on whether an 
oscilloscope sweep of 50 ns or 100 ns is used), and a linear interpolation near the 
minimum is used to obtain t precisely. 

Another method uses the symmetry of the received pulse by finding the point t^ 
that divides the area under RP(t) into two equal halves. 


The correlation method does not give a center of the received pulse, but instead 
determines the time T that the received pulse must be shifted back toward the trans- 
mitted pulse to achieve the best correlation between the two pulses. Since 
T cl - T rc - T tc and it is this difference in pulse centers, rather than their individual 
values, that is used in determining the range correction A, the method is applicable 
to the purpose at hand. The value of T CL is determined as that value that maximizes 
the cross-correlation function, i. e., the value that gives 


J 


TP(t) • RP(t + T cl ) - max 


( 7 ) 


In practice, the integration is replaced by summation, and T 
two discrete values. 


is interpolated between 


Once T.^ and T^,^, or their difference has been determined, it remains to deter- 
mine either the threshold values T,^ and T^, or T^,^. - T^,. Since the oscilloscope 
sweeps are started by trigger pulses from the counter’s START and STOP circuits, 
T^rp - T-™ is constant to within the jitter of the electronic circuits involved. If this 
constant has been included in the calibration of the system’s electronic delay, the 
range correction becomes T pr - T™^. 


If the jitter is appreciable, it can be eliminated from the range correction by 
determining T^,^, and T R( j, from the corresponding threshold voltages V r p and V R . 
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Although it is difficult to measure these voltages directly because they are somewhat 
dependent on the slope of the pulse's leading edge, they can be obtained experimentally 
by ranging to a reflector in a fixed position. The method, which works best when 
there are large pulse-to-pulse variations in the shape of the leading edge, is the 
following: a series of threshold voltages are assumed, and the corresponding threshold 
times are determined from the photographs for each pulse in the series. Then the 
sample mean and the sample variance of the threshold time are computed for each 
assumed voltage. The threshold corresponding to the smallest sample variance cor- 
responds to the actual threshold. Without jitter, the threshold should occur at a fixed 
position on the x axis of the oscilloscope and the variance should be zero. Since there 
is jitter, the sample variance will not be zero at the actual threshold voltage, but equal 
to the rms value of the jitter between the threshold activation of the counter and the 
start of the sweep. 

When the values of V T and V R are precisely known, T^ and T RT can be deter- 
mined from the oscilloscope photographs, and the jitter between counter and oscilloscope 
does not enter into the measurement of the range correction. There is a method that 
can be used to determine V T and V R without assuming that T TT and T RT have fixed 
positions on the oscilloscope photographs. This method uses range measurements 
made to a reflector in a fixed position. A series of threshold voltages are again 

assumed for the transmitted and received pulses. For the measurements correspond- 

* 

ing to each pair of assumed thresholds, the range correction A is calculated and added 
to the counter reading. Then, the sample mean and the sample variance are calculated 
for all the fixed- reflector measurements. Because this range is fixed, the variance 
is a minimum when the assumed thresholds are equal to the actual thresholds. This 
minimum variance is a measure of the range-correction error for the received signal 
strength {in photoelectrons) that was used for the experiment. 

4. RESULTS 

Tables 1 and 2 summarize initial results obtained from range measurements made 
in accordance with the methods described in previous sections. These results are 

* 

In practice, the transmitter threshold voltage is determined with sufficient accuracy 
by the previously described procedure, and only variations in the receiver threshold 
are assumed. 
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illustrative of the procedures involved- They show that the corrections made from the 
oscilloscope photographs do improve the accuracy of the counter readings. The results 
are consistent with those expected from the equipment in current use at the SAO observ- 
ing stations, but they are not representative of the ultimate accuracy to be expected 
eventually with improved equipment and procedures. It can be seen from Table 1 that 
or the jitter between counter and oscilloscope, is about 1 ns, with most of this error 
associated with the received pulse. Although the counter's threshold discrimination 
and trigger-generating circuits are identical in both START and STOP channels, they 
respond differently to the smooth leading edge of the transmitted pulse than they do to 
the irregular leading edge of the received pulse. Hence, it appears reasonable that 
the jitter associated with the received pulse is larger than that associated with the 
transmitted pulse. 


The range error applies to range corrections obtained by setting A = T RC - T T £. 

These corrections are based on the assumption that T™™ - T„_, = const. Hence, ff . 

1 1 Ki 1 

includes the effect of the synchronization jitter. Table 1 shows that cr^ is about 1. 5 ns. 
The quantity a % applies to measurements made by using the correction 

A= ^ T RC " T RT* “ ( T TC " T TT^ V T and V R selected t0 minmuze a 2 - The a g 
values in Table 1 are about 1 ns, a value that appears to be reasonable for a return 

pulse of 500 electrons. For a return of this strength, we would expect an error of 

1 /2 

8. 5/500 ' “ = 0. 4 ns in determining the center of the pulse and an additional error of 
1/2 * 

1/12 = 0. 3 ns for a counter with a uniform probability density over its 1-ns resolu- 

tion interval. 


Table 2 shows that the residuals of three satellite transits were reduced by apply- 
ing the corrections obtained from oscilloscope photographs. Range corrections for 
several other satellite transits were also computed, but improvements in the short- 
arc residuals could not be realized because in some cases the variance of the uncorrected 
range residuals was considerably larger than that of the range corrections. In one case, 
the average number of electrons per pulse was too small to make an effective determina- 
tion of the center of the received pulse. 

The satellite BE-C shows the greater decrease in residual (0. 7 to 0.4 m), a result 
that might be expected because the retroreflectors on this satellite are arranged in such 
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a way that they produce less of a random variation in range with changes in the satel- 
lite’s aspect angle than do those on Geos 1. For the measured pulse width of the 

BE-C returns, the corresponding a of a normal pulse is about 8. 5 ns; hence, the error 

1/2 

in determining the center of a 38-electron pulse is 8. 5/38 = 1. 4 ns (0. 2 m). This 

error combined with the 1. 7-ns synchronization jitter in Peru gives 
(1.4 2 + 1. 7 2 ) 1//2 = 2. 2 ns (0. 33 m), a value comparable to the 0. 39 m shown in 
Table 2. 

In computing the short-arc residuals shown in Table 2, a 4- to 5-day orbit was 
first obtained by using laser data from all stations. Then, best -fitting short arcs 
for the given pass were obtained by varying inclination, eccentricity, mean anomaly, 
and mean motion. All range observations used in the orbit computations were corrected 
for the nonrandom displacements between pulse center and the center of mass of the 
satellite. These corrections, which depend on the angle between the laser beam and 
the axis of the satellite's retroreflector, were computed by Arnold (1972). 
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Table 1. Summary of range measurements to fixed reflectors. 


Station 

Oscilloscope 

sweep 

duration 

(ns) 

Number 
of pulses 

Average 
number of 
electrons 
in received 
pulse 

Measured pulse 
durations 

Synchronization jitter 


Range error 

Transmitted 

pulse 

(ns) 

Received 

pulse 

(ns) 

Transmitted 

(ns) 

Received 

(ns) 

a 

c 

(ns) 

a i 

(ns) 

(ns) 

Brazil 

50 

17 

560 

22 

26 

0. 1 

1.3 

2. 1 

2. 1 (CT) 

0.7 (CL) 

Peru 

50 

18 

463 

2G 

32 

0.2 

1.3 

1.7 

1. 5 (CT) 

1. 1 (CT) 

Arizona 

100 

20 

416 

21 

26 

0.3 

0. 9 

1.7 

1.3 (S) 

1.3 (CL) 


The a’s are sample standard deviations of the series of pulses: cr represents the uncorrected counter 
readings, and o\ and <7 represent corrected counter readings. For o the counter and oscilloscope were 
assumed to be in exact synchronism. For the receiver threshold was determined by using the fact that 
the range is constant. The cr and <7 values given above are the lowest for the four methods of pulse-center 
determination: centroid (CT), symmetry (S), half-area (A), or correlation (CL). 



Table 2. Summary of satellite -range corrections. Method used for pulse-center determination: 
centroid (CT), symmetry (S). 


Date 


Oscilloscope 

sweep 

Station duration 

location Satellite (ns) 


Number of 
observations 


Average 
number of 
electrons 
in received 
pulse 


Measured pulse 
widths 

Transmitted Received 
(ns) (ns) 


Short-arc residuals 

Uncorrected Corrected 

observations observations 
(m) (m) 


7 June 1972 

Brazil 

BE-C 

50 

9 

17 June 1972 

Brazil 

Geos 1 

50 

35 

17 March 1973 

Peru 

BE-C 

100 

19 


38 

24 

22 

0.73 

0.45 (CT) 

85 

20 

21 

0.71 

0.57 (S) 

276 

23 

19 

0.74 

0. 39 (CT) 
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Figure 1. Block diagram for pulse detection. 
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Figure 2. Diagram for the range correction A, where T R = T c + A and A= (T RC - T RT ) - (T TC - T TT ). 


Digitized representation of the photographed pulses: * = transmitted pulse; x = received pulse, 
duration = 100 ns. 


Figure 3. 


Sweep 



